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Abstract 

Rapidly rotating, stably stratified three-dimensional inviscid flows conserve both energy and 
potential enstrophy. We show that in such flows, the forward cascade of potential enstrophy 
imposes anisotropic constraints on the wavenumber distribution of kinetic and potential energy. 
The horizontal kinetic energy is suppressed in the large, nearly horizontal wave modes, and should 
decay with the horizontal wavenumber as /c^ 3 . The potential energy is suppressed in the large, 
nearly vertical wave modes, and should decay with the vertical wavenumber as A:" 3 . These results 
augment the only other exact prediction for the scaling of energy spectra due to constraints by 
potential enstrophy obtained by Charney (J. Atmos. Sci. 28, 1087 (1971)), who showed that in 
the quasi-geostrophic approximation for rotating stratified flows, the energy spectra must scale 
isotropically with total wavenumber as fc -3 . We test our predicted scaling estimates using resolved 
numerical simulations of the Boussinesq equations in the relevant parameter regimes, and find 
reasonable agreement. 
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Classical quasi-geostrophic (QG) flow is a useful approximation for rapidly rotating, 
strongly stratified flows Q, Q]- In this approximation, the zeroth-order expansion of the 
velocity in the rotation and stratification parameters is geostrophic, meaning that the Cori- 
olis force is balanced by the pressure gradient force. Furthermore the linear plane- waves 
called inertia-gravity waves are eliminated to the lowest order. This leads to a simplification 
of the dynamics which is described entirely by the evolution of potential vorticity q qg j^]: 

+ U 0h ■ Vq qg = 0, (1) 

d6 

where q qg = f— - Nu 3 , (2) 

where u 0h is the leading order (horizontal) geostrophically balanced velocity, 9 is the den- 
sity fluctuation scaled to have the same dimensions as velocity, 003 = (Vh x Uh) ■ z is the 
z-component of the vorticity, / and N are the Coriolis and Brunt- Vaisala (buoyancy) fre- 
quencies respectively of a system which is rotating stratified in the z-direction. In QG flow, 
the vertical velocity w appears as a correction to leading order geostrophic balance. 

In 1971, Charney [3|] showed that the global conservation of total energy — \ J (\u\ 2 + 
9 2 )dx, and potential enstrophy Q qg — ~ J \q qg \ 2 dx by inviscid three-dimensional (3d) QG 
flow is analogous to conservation of energy and enstrophy in non-rotating two-dimensional 
turbulence (see for example 4J for further discussion on the assumptions and details of 
Charney's work). Following the classical theory of 2d turbulence {5, 6], Charney predicted an 
inverse (upscale) cascade of energy with corresponding large-scale energy spectrum Ef{k) oc 
£ 2 / 3 /c -5 / 3 , and a forward (downscale) cascade of potential enstrophy with corresponding 
small-scale energy spectrum Ex{k) oc £g 3 fc -3 , where k is the wavenumber and e and Eq are 
the transfer rates of energy and potential enstrophy, respectively. 

In theory, as the rotation and stratification of a 3d fluid become infinitely strong, the 
inertia-gravity waves are eliminated to lowest order, giving leading order QG flow satisfying 
( TT|2|) . In practice, for very large but finite rotation and stratification, the inertia-gravity 
waves strongly influence the small-scale dynamics leading to Ex{k) ~ /c~ 7 where 1 < 7 < 2 



in the high wavenumbers. The underlying leading order QG scaling of k~ 3 predicted 
by Charney can then only be extracted by separating the QG (or geostrophic) modes from 
the wave (or ageostrophic) modes by either suitably projecting the full solutions onto the 
QG modes 0,0] or by filtering out the ageostrophic inertia-gravity waves [l|. 

In the present work we consider rotating stratified turbulence retaining both the leading 



order QG as well as all sub-leading contributions from inertia-gravity waves and other non- 
linear waves. We show that in this parameter regime, as in classical QG, potential enstrophy 
plays a significant role in constraining the energy. For a wavevector k = k x x + k y y + k z z 
the horizontal component is kh = (k 2 + k 2 ) 1 ^ 2 and the vertical component is k z . We show 
that the potential enstrophy dominates over potential energy in the large, nearly vertical 
modes (k z /kh 1), resulting in a potential energy spectrum scaling of k~ 3 for large k z . 
And potential enstrophy also dominates over horizontal kinetic energy in the large, nearly 
horizontal modes (kh/k z ^> 1), resulting in a horizontal kinetic energy spectrum scaling of 
kfr 3 for large kh- These are the first scaling estimates for the spectra of rapidly rotating and 
stably stratified flows away from pure QG, obtained solely using the relationship between 
potential enstrophy and energy as a function of scale. Such scaling laws, apart from being 
benchmarks, are potentially useful in parameterizing the turbulent small scales in large sim- 
ulations of rotating stratified flows thus reducing the computational expense of explicitly 
resolving the small scales. 

We begin with the Boussinesq equations for rotating, stably stratified and incompressible 
flow given by flo| : 

—u + fz x u + Vp + N9z = uV 2 u + T 

—e-Nw = kV 2 9 (3) 
Dt 

Vu = 0, 

D d 

where — = — +u- V, u is the velocity, w is its vertical component, p is the effective pressure 
Dt dt 

and T is an external input or force. The total density is pr{x) = po — bz + p(x), where po 
is the constant background, b is also constant and larger than zero for stable stratification 
in the vertical z-coordinate, p is the density fluctuation such that \p\ \bz\ po an d 
9 = p(g /bpo) 1 / 2 has the dimensions of velocity. The Coriolis parameter / = 2Q where Q is 
the constant rotation rate about the z-axis, the Brunt- Vaisala frequency iV = {gb/po) 1 ^ 2 , 
v = p/po is the kinematic viscosity and k is the mass diffusivity coefficient. We assume 
periodic or infinite boundary conditions. The relevant non-dimensional parameters for this 
system are the Rossby number Ro = f n i/f and the Froude number Fr = f n i/N, where 
fni — ( e /^/) X ^ 3 is the non-linear frequency given input rate of energy e/ Thus Ro and Fr 
are the ratios of rotation and stratification timescales respectively to the nonlinear timescale. 
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The Boussinesq equations conserve the following quantities for T = v = k 

D f 

total energy E T = E + P, — / E T dx = 0, 

Dq 



potential vorticity q = (u> a ■ V pxj , = 



potential enstrophy Q = -q 2 , — p = — f Q dx = 0. 

E = \\u\ 2 is the kinetic energy, P = ^9 2 is the potential energy of the density fluctuations. 
The absolute vorticity uj a = uj + fz and the relative (or local) vorticity uj = V x u. Potential 
vorticity may be written in terms of 9 as 

an 

q = fN + u-V6 + f—-Nuj 3 . (4) 

oz 

The constant part fN does not participate in the dynamics and we will therefore neglect it 
from now on. The linear part of is precisely q qg of (j2J). In what follows we will assume 
that v — > and k — > such that Prandtl number Pr — v/k — 1, and the force JF is 
confined to the lowest modes. Thus we assume a conventional 'inertial-range' of turbulent 
scales wherein the transfer of conserved quantities dominates over both their dissipation and 
forcing. 

A 8 ^ ^ and Fr 0, the potentia. vorticity , approaches fa Q Q. Th,s is easi.v 
observed by considering the non-dimensional form of (j3J) namely q = uj ■ W + i?o _1 || — 
Fr _1 ti;3, and letting i?o and Fr tend to zero together. In fourier representation: 

q{k) ~ + iA^fe ft x n h = /A;^ + ziVA;^w ft (5) 

where 7 denotes fourier coefficients, the total wavevector k = kh+k z z, the horizontal wavevec- 
tor component has length kh = (k 2 + k 2 ) 1 ^ 2 , the vertical wavenumber is k z and Uh is the 
horizontal velocity vector with magnitude Uh = (u 2 + u 2 ) 1 ^ 2 . We assume that the vertical 
velocity w = u z ~ in the lowest order (classical QG) thus obtaining the last equality of 
Eq. ©. 

We take both N and / to be very large, and N/ f = 1 so that Ro = Fr. This approaches 
the special case Ro — > and Fr — > while i?o = Fr which was shown rigorously to be 
leading order QG in [h]]. For k > kf, scales smaller than the forcing scale, we consider two 
cases: 

k 

1 ) Case ^> 1 . These are the more vertical wavenumber modes corresponding loosely to 
flat 'pancake' scales in physical space. Eq. (j3J) reduces to q ~ fk z 9, yielding the following 



relation between potential enstrophy and potential energy distribution in spectral space, 

Q(k h ,k z ) = ^\q\ 2 = f 2 k 2 z P(k h ,k z ), 

which upon integration over some high vertical wavenumber interval leads to the following 
constraint: 

/•OO POO 

Q{k h , k z )dk z > ] 2 k\ I P(k h , k z ) dk z , 



where the potential energy spectrum P(kh, k z ) = \\d\ 2 . Thus, for sufficiently high wavenum- 
bers k z — > oo, the potential enstrophy Q forms the dominant forward cascade and, in order 
to remain finite, suppresses the potential energy P in this regime. The dimensional argu- 
ment following p, |6|] assumes that in this wavenumber limit, the potential energy spectrum 
must depend on the potential enstrophy flux rate eq and the vertical wavenumber k z , so 
that: 

P(k h , k z ) ~ efkf. (6) 

kh 

2) Case -f- > 1 . These are the wide flat wavenumber modes corresponding to the tall 

k z 

columnar scales in physical space. In this limit Eq. (jSD reduces to q = iNkhith- Following 
the same arguments as for potential energy above, we obtain that the potential enstrophy 
dominates the forward cascade in the regime kh/k z 3> 1, resulting in suppression of horizontal 
kinetic energy resulting in the following scaling estimate: 

E h (k h ,k,)~e% B k ?. (7) 

where the horizontal kinetic energy Eh{kh, k z ) = ^\uh\ 2 - 

Our new scaling predictions (jSJ) and © describe the spectral statistics of the full flow for 
small but finite Ro and Fr, not just the leading order QG part of the dynamics. The main 
purpose of this paper is to show that the potential enstrophy imposes predictable constraints 
on the energy in rapidly rotating stratified flow even when the flow is not strictly QG. In 
the next sections we seeking to numerically verify our predictions for the energy spectra in 
the two limiting regimes in wavenumber. 

We simulate the Boussinesq equations taking N = f very large. We use a pseudo- 
spectral code in a periodic cube of side L — 1, generating wavenumbers which are integer 
multiples of 2n. A fourth-order Runge-Kutta time integration is used and the inertia-gravity 
wave frequencies are resolved in our explicit scheme. Since we are interested in the small 
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# n k f N/f Ro Fr e f 

1 256 4 1 0.029 0.029 0.62 

2 256 4 1 0.014 0.014 0.61 

3 256 4 1 0.0072 0.0072 0.60 

4 512 4 1 0.0072 0.0072 0.60 

TABLE I: Parameters of Boussinesq calculations: n - number of grid points to a side; kf - forcing 
wavenumber; N - Brunt- Vaisala frequency, Ro - Rossby number; Fr - Froude number; ef - rate 
of input of kinetic energy. 

scales (high wavenumbers) , we use a low- wavenumber stochastic forcing centered at kf — 4. 
In order to extend the inertial range in the high wavenumbers, the viscous dissipation is 
modeled using a hyperviscous term (— l) p+1 h>(V 2 ) p u, where p = 8 in place of the normal 
laplacian viscosity term z/V 2 it. The hyperviscosity coefficient v is dynamically chosen based 
on the energy in the highest mode for both momentum and mass diffusion following Q], 

1 /2 

v{t) = 2.5 ^ £, (^ m ' t ) ^ k 2 m 2p where k m is the highest available wavenumber and E(k m ,t) is 
the kinetic energy in that wavenumber. An analogous scheme is used for the diffusion term 
in Eq. for the evolution of 9. The parameters of a few of our runs are given in Table [B 
We report the results from data #4 for which Ro and Fr are the smallest and the resolution 
is the best. 
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FIG. 1: Total potential enstrophy Q and its linear part Q qg for data #4. The two are indistin- 
guishable, indicating purely quadratic potential enstrophy. Time t is in dimensional units. 

Figure [1] shows the evolution of the total potential enstrophy Q from (j3]), and its linear, 
quasi-geostrophic piece Q qg . The former is indistinguishable from the latter indicating that 
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the nonlinear part of the potential vorticity oj ■ V6 is negligible and thus the potential 
enstrophy is quadratic. The system is QG in the leading order, or near-QG in the sense 
described above. The mean potential enstrophy has reached a nearly steady value in the 
time range 2.5 < t < 5.2, which corresponds to between 5 and 11 non-linear time cycles, or 
5000 to 10400 rotation (stratification) cycles. 



FIG. 2: Log-log plot of spherical shell averaged potential and kinetic energy spectra for data #4 
at time t = 5.2. The high wavenumber scaling is /c _1 indicating that in this representation the 
energy is dominated by waves. Inset: Same spectra averaged over time 3 < t < 5.2. 

Figure [2] shows the shell-averaged kinetic and potential energy spectra for our simulation, 
computed as follows: 



where A; — 0.5 < k' < A; + 0.5 thus including all wavenumbers in the spherical shell of average 
radius k. The scaling of both E(k) and P(k) is k~ l for k 3> kf which indicates that by this 
measure the high wavenumbers are still dominated by waves [l2]. The inset shows the 
average of the spectra over the time period 3 < t < 5.2 over which the potential enstrophy 
is constant as seen in Fig. [TJ The time-averaged spectra also show a scaling very close to 
k~ x indicating that the small scales (k > kf) have achieved close to a statistically steady 
state. 

The potential energy and horizontal kinetic energy spectra as functions of kh and k z were 
computed as double sums according to: 




k 



E(k) = -J2Hk')\", P(k) = -J2\0(k')\ 2 



k' k' 
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FIG. 3: Log-log plot of potential energy density P(kh,k z ) vs. k z for data #4 averaged over time 
3 < t < 5.2. Each curve is the spectrum for a different fixed value of kh- For 10 < kh < 50 and 
k z S> A;/, the scaling ranges between k~ 4 and /c~ 3 . Inset: Same spectra versus k z /kh shows that 
the 'turnover' to the inertial range scaling for all the curves emerges only when k z /kh > 1. 

where k z [k h ] — 0.5 < k' z [k' h ] < k z [k h ] + 0.5. Figure [3] shows P(kh,k z ) as a function of k z 
for various values of kh- For 10 < kh < 50 and kf < k z < 100, the scaling for P(kh,k z ) 
ranges between fc~ 4 and /c~ 3 indicating stronger suppression of potential energy than the 
dimensional prediction of Eq. ([6]). As kh increases P(kh, k z ) also persists more strongly into 
the high k z . Conversely, for a fixed small k z < kf, the smaller kh spectra have more energy, 
indicative of a growth of potential energy as k z <C kf for small kh- The inset of Fig. [3J 
shows the same spectra versus k z /kh which shows clearly that the inertial range scaling 
for each kh emerges only when k z /kh > 1, the predicted range for Eq. ([6]). Overall the 
constraints on potential energy due to potential enstrophy are thus highly anisotropic in 
scale and consistent with our prediction. 

Figure H] shows Ehkh,k z as a function of kh for various values of k z . For 10 < k z < 50 
and 10 < k h < 100, the horizontal energy spectrum E h (kh,k z ) scales between k^ 4 and fc^ 3 
consistent with the suppression of horizontal kinetic energy by potential enstrophy in these 
modes. As kh grows, the horizontal kinetic energy persists more strongly into the high kh- 
Conversely, for a fixed small kh, the smaller k z have more energy, indicating a growth of 
energy upscale in k z . The inset of Fig. H] shows again that the inertial range scaling predicted 
arises only in the anisotropic regime kh/k z > 1 consistent with prediction. 

In conclusion, we have deduced separate scaling laws for horizontal kinetic and potential 
energy spectra, both of which are constrained in the small scales by potential enstrophy in 
rapidly rotating stably stratified flows. Potential enstrophy suppresses the potential energy 
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FIG. 4: Log-log plot of horizontal kinetic energy density Eh(kh,k z ) vs. kh for data #4 averaged 
over time 3 < t < 5.2. Each curve is a different value of k z . For 10 < k z < 50 there the scaling 
ranges between k~ 4 and k~ 3 . Inset: Same spectra vs. kh/k z showing that the inertial range 
emerges as kh/k z > 1 as predicted. 

in the large, nearly vertical modes, and also suppresses horizontal kinetic energy in the 
large, nearly horizontal modes; the resulting energy densities in these modes scales as k~ 3 
and /c^ 3 respectively. Our test simulations data show even steeper scaling of the spectra 
than predicted (greater suppression due to potential enstrophy). The numerical calculations 
used to verify our predictions are, at 512 grid-points to a side, the highest resolution unit 
aspect-ratio simulations of the Boussinesq equations performed to date. Higher resolution 
may well show closer agreement with our theoretical prediction; we have already observed a 
tendency toward our predicted —3 exponent when going from 256 3 to 512 3 in grid-resolution. 
Most importantly, the scalings predicted and observed are very different from the isotropic 
/c~ 7 (1 < 7 < 2) scaling expected, for example, in the wave-dominated shell-averaged spectra 
for near-QG flows (see Fig. [T]). Our only assumption is that rotation and stratification are 
strong enough that the potential vorticity becomes linear, and hence the potential enstrophy 
quadratic. We do not invoke additional asymptotics nor do we need to limit ourselves only 
leading order modes. In future work we will extend our analysis to the case of N/f ^ 1, 
that is, the strength of rotation and stratification are large but unequal. The possibilities 
for generalized quasi-geostrophic flow [13| in which aspect ratio is an additional parameter, 
are also promising areas for further research. 
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